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I  -  ABSTRACT 

We  report  on  progress  made  on  USAFOSR  Grant  No.  AF0SR-80-0083.  Substantial 
results  on  the  use  of  mixed  finite  element  methods  for  partial  differential 
eigenvalue  problems  and  for  viscous  flow  problems  have  been  achieved.  Both 
analytical  and  computational  studies  have  been  carried  out.  The  analytical 
work  is  concerned  with  error  estimates  for  the  finite  element  approximations, 
while  the  computational  efforts  consist  of  implementing  these  finite  element 
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II  -  DESCRIPTION  OF  WORK 

In  our  proposal,  we  proposed  to  consider  the  approximate  solution, 
via  mixed  finite  element  methods,  of  two  physical  problems.  The  first  was 
the  energy  stability  of  the  incompressible  Navier- Stokes  equation.  This  problem 
has  the  following  variational  formulation:  seek  (tj,  <p,  v*)  €  V  *  S  *  H1 
such  that 


(£  •  D  •  w  -  (j,  divw)  =  -v*  •  vw 

(2.1) 

* 

ip  div  =  0 


for  all  (w,  *)  €  V  x  S  .  Here  V  and  S  are  appropriately  chosen  Hilbert 
spaces  and  D  is  the  deformation  tensor  of  a  given  flow.  It  is  easily  shown 
[1]  that  if  the  largest  eigenvalue  v  of  the  problem  (2.1)  is  smaller  than  the 
kinematic  viscosity  v  of  the  given  flow,  then  the  given  flow  is  stable  in 
the  energy  sense.  We  note  that  (2.1)  is  a  linear  eigenvalue  problem  which 
determines  stability  regions  for  arbitrary  perturbations,  i.e.  not  infinitesimal. 
The  problem  (2.1)  is,  however,  not  easily  solvable.  We  therefore  wish  to 
solve  the  problem  approximately.  To  this  end  we  choose  finite  dimensional 
subspaces  /cV  and  Sh  c  S  and  require  that  (2.1 )  hold  for  all  (w,  ^) 

€  V*1  x  .  This  leads  to  a  linear  generalized  algebraic  eigenvalue  problem 
of  the  form 


(2.2) 


We  choose  the  spaces  V*1  and  Sh  to  be  finite  element  spaces. 
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Under  certain  conditions  on  these  spaces,  conditions  which  are  necessary  for 
the  stability  of  the  approximate  problem,  we  have  been  able  to  prove  that  the 
approximate  eigenvalues  converge  to  the  exact  eigenvalues  v*  at  an 
optimal  rate.  In  addition,  two  computer  codes  have  been  written  which  use 
different  finite  element  spaces  which  satisfy  the  above  mentioned  conditions. 
The  first  uses  a  tri angulation  of  the  type  illustrated  in  the  figure  on  the 
left  wherein  the  vector  is  approximated  by  piecewise  linear  polynomials  in 
each  triangle  and  the 


Figure  1 

scalar  $  is  approximated  by  piecewise  constant  functions  in  each  triangle. 
The  second  code  uses  the  element  on  the  right  which  again  uses  piecewise 
linear  vectors  on  triangles,  but  now  uses  scalars  which  are  constant  through- 

L. 

out  the  quadrilateral .  Our  theory  predicts  that  v  converges  to  v* 
quadratically  in  the  grid  spacing  and,  of  course,  this  is  verified  by  our 
computational  results.  The  codes  developed  take  advantage  of  the  structure 
and  sparsity  of  the  algebraic  problem  (2.2). 

The  second  physical  problem  considered  is  that  of  acoustic  eigenvalue 
problems  such  as  those  described  by  the  equation 


Ad>  =  X4>  in  n 


(2.3) 
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Again  this  eigenvalue  problem  can  be  given  the  variational  characterization: 


seek  (u,  <j>.  A)  e  V  *  S  *  F  such  that 

il)  div  u  *  J 

■’fi  h 

(2.4) 

{ div  v^  +  v^  •  u_)  =  0 

for  all  (v,  ♦)  f  V  k  S  ,  where  again  V  and  S  are  suitable  Hilbert  spaces. 
Approximations  are  defined  as  before.  In  this  case  we  have  proved  optimal 
error  estimates  for  the  eigenvalues  X  and  ei genfunction  pairs  (^,  4>)  and 
have  developed  codes  to  compute  these  eigenvalues  and  eigenfunction,  again 
using  the  two  elements  described  in  the  figure  above. 

The  variational  problems  (2.1 )  and  (2.4)  are  special  cases  of  the  abstract 
mathematical  problem:  seek  (u,  4,,  x)  €  V'1  *  S1  x  £  such  that 


a(u,  y_)  +  b^v,  41)  =  Ac(^,  v) 
b2(iu  i|>)  =  xd(ip,  X) 


(2.5) 


for  all  (vs  *)  {  k  ,  Here  a,  b^ ,  b2,  c  and  d  are  sesqui linear  forms 
defined  on  the  appropriate  Hilbert  spaces.  We  have  been  able  to  prove  optimal 
error  estimates  for  eigenvalue  problems  of  the  type  (2.5)  for  a  variety  of 
combinations  of  forms  a,  b^ ,  etc.  under  reasonable  hypothesis  on  these  forms. 
For  example,  the  cases  b^  =  b2  with  c  or  d  =  0  ,  b^  =  b^  with  both  c 
and  d  f  0  ,  and  b^  /  b^  ,  all  with  a(iu  vj  coercive  or  weakly  coercive  have 
been  successfully  analyzed  and  illustrative  computer  programs  have  been 
implemented  in  each  case.  The  various  cases  which  have  considered  describe 
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a  variety  of  physical  examples  includinn  some  generalized  non-self-  adjoint 
acoustic  eigenvalue  problems,  eigenvalue  problems  emanating  from  transmission 
line  theory,  etc. 

Most  of  the  work  described  above  concerning  the  use  of  mixed  finite 
element  methods  for  eigenvalue  problems  constitutes  the  Ph.D.  dissertation  of 
Ms.  Janet  Peterson,  whose  studies  have  been  supported  by  the  AFOSR  grant. 

This  work  was  reported  on  at  the  SIAM  Fall  1980  meeting  in  Houston  and  will 
be  the  subject  of  a  series  of  papers  now  in  preparation. 

In  our  proposal,  we  also  proposed  to  develop  mixed  finite  element 
discretizations  which  satisfy  certain  conditions  (described  in  our  previous 
proposal)  which  insure  the  stability  and  optimal  accuracy  of  the  discretizations. 
We  have  carried  out  this  study,  concentrating  on  second  order  L,, 
approximations  to  the  velocity  vector  in  viscous  flow  problems.  So  far  our 
study  has  considered  only  the  linear  Stokes  or  Oseen  equations.  The  elements 
considered  are  those  of  Figures  1  and  2. 


In  Figure  2,  on  the  left,  the  velocities  and  pressures  are  approximated  on 
each  quadrilateral,  by  bilinear  and  constant  function,  respectively,  while 
for  the  element  on  the  right,  the  velocities  and  pressures  are  approximated, 
on  each  triangle,  by  piecewise  quadratic  and  piecewise  constant  functions. 
This  last  element  appears  often  in  the  literature  [2].  For  all  four  elements 


the  velocity  approximations  are  second  order  (in  the  grid  spacing)  accurate  in 
l_2  while  the  pressure  approximations  are  first  order  accurate.  However,  the 
elements  pictured  on  the  right  of  Figure  1  and  on  the  left  of  Figure  2 
achieve  this  rate  of  accuracy  with  considerably  fewer  unknowns  than  the  other 
two  elements.  Computer  codes  implementing  all  four  elements  have  been  written 
and  are  presently  being  used  to  generate  information  which  will  be  used  to 
compare  the  elements. 

This  work  was  also  reported  on  at  the  SIAM  Fall  meeting  and  in  various 
papers  (1,  2,  3,  4  of  Section  3). 

[1 ]  Serrin,  J.:  "Mathematical  principles  of  classical  fluid  mechanics", 
Handbuch  der  Physick,  8(1959),  Springer. 

[2]  Girault,  V.  and  P.  Raviart:  Finite  Element  Approximations  of  the 
Navier-Stokes  Equation,  (1979),  Springer. 
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